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We present an exact solution of the Worm-Like Chain (WLC) model for semi-flexible polymers valid over the entire range of 
polymer lengths. Our results are in excellent agreement with recent computer simulations and reproduces important qualitatively 
interesting features observed in simulations of polymers of intermediate lengths. We also make a number of predictions that 
can be tested in a variety of concrete experimental realizations. The expected level of finite size fluctuations in force-extension 
curves is also estimated. This study is relevant to mechanical properties of biological molecules. 



PACS numbers: 87.45-k,05.40+j,36.20-r 

Many biologically important molecules, like DNA and 
Actin, are semi- flexible polymers jy. In recent years 
there have been experiments j^] which pull and stretch 
single molecules to measure elastic properties. For in- 
stance, one can study the "equation of state" of a semi- 
flexible polymer by measuring its extension Q as a func- 
tion of applied force. Alternatively one can tag the ends 
with fluorescent dye Q and determine the distribution 
of end-to-end distances. Such studies reveal a wealth 
of information about the mechanical properties of semi- 
flexible polymers, which is of clear biological importance. 

Till a few years ago, studies of polymer molecules such 
as DNA were limited to samples containing large num- 
bers of molecules ||. This made it hard to probe the 
elastic properties of individual DNA which are of vi- 
tal importance to biological processes such as protein- 
induced DNA bending ||||. It is only quite recently, 
due to advances in technology that single molecule stud- 
ies became feasible. In order to correctly interpret single 
molecule experiments which are now being performed, 
a good theoretical understanding of semi-flexible poly- 
mers is essential. Quite apart from the biological interest, 
semi- flexible polymers are of interest to physicists j^U ■ 
This paper is devoted to understanding the equilibrium 
statistical mechanics of single semi-flexible polymers. 

Statistical mechanics of a single polymer molecule is 
dominated by fluctuations because it is a system of fi- 
nite size. It is only in the thermodynamic limit of ex- 
tremely long polymers that these fluctuations about the 
mean die out. Due to the dominance of fluctuations, the 
experimentally measured mean values for a semi-flexible 
polymer crucially depend on the precise choice of the en- 
semble. For instance, one gets qualitatively distinct fea- 
tures in force-extension curves depending on whether the 
force or the extension is held constant in an experimental 
setup |@. 

The most popular theoretical model for understanding 
semi-flexible polymers is the Worm-Like Chain (WLC) 
jTlJ, which ignores self-avoidance and models the poly- 
mer as a framed space curve of fixed total length L with 
an energy cost for bending and twisting. In order to in- 



terpret the experimental data, it would be useful to have 
a clear and complete understanding of the predictions of 
the WLC model. Such an understanding would reveal 
the strengths and deficiencies of the model in describing 
real polymers and could be used to improve the model. 
There do exist partial results jl2| on the statistical me- 
chanics of the WLC model: some theoretical studies 0] 
investigate the flexible limit of very long polymers (long 
compared to the persistence length;!/ >> L p ). However 
experimental interest is not confined to very long poly- 
mers. For example, experiments on Actin |4j deal with 
polymers of length L = 30/im, which is only about twice 
the measured persistence length of L p = 16.7/ito. There 
is also theoretical work |]l2|| on extremely short polymers. 
There is a clear gap in the present understanding of the 
WLC model for polymers of intermediate length. Our 
purpose in this paper is to fill this gap. We present a 
solution of the WLC model and describe the equilibrium 
elastic properties expected from the model. Our solution 
is exact in the sense that the elastic properties can be 
determined to any desired accuracy. The main results of 
this paper are contained in the figures which show the 
force-extension relation and end-to-end distance distri- 
butions predicted by the WLC model. These predictions 
agree well with two independent computer simulations 
g|§ (Figs.3 and 4). 

WLC model with pure bend: A configuration C of the 
polymer is described by a space curve x(s), with s the 
arc-length parameter (0 < s < L) ranging from to L, 
the contour length of the polymer. The tangent vector 
t = dx I ds to the curve is a unit vector 

t.i=i (i) 

and the curvature of the polymer is given by k — \dt/ds\. 
We will suppose that one end of the polymer is tethered 
to the origin (x(0) — 0) and the other end x(L) = r is 
tagged. As the polymer configuration changes with ther- 
mal agitation, the location r of its tagged end fluctuates. 
The quantity we wish to compute is Q(r), which is the 
probability distribution for the location r of the tagged 
end 0. If the tagged end is pulled from r to f+dr, Q(f) 
changes and consequently, the free energy. This implies 
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that a force is needed to stretch the polymer. Thus Q(r) 
is directly related to the force-extension relation of the 
polymer. To compute Q(r) we need to sum over all poly- 
mer configurations C which end at r, with a Boltzmann 
weight: Z = Sc cxp(— £[C]//csT), where the energy £ 

associated with a configuration C is £{C) = \A J Q L dsK 2 
and A is the bending modulus. This is a standard count- 
ing problem in statistical mechanics and can be natu- 
rally addressed in the language of path integration [fl3|| . 
However, not much progress has been made because of 
the difficulty (ljjl^l presented by the inextensibility con- 
straint ([!]). The key to circumventing this difficulty is 
to consider [§];[l6) Brownian motion in the space of tan- 
gent vectors (t) rather than (as is customary for flexible 
polymers) position vectors x. The tangent vectors form 
a unit sphere (See Eq. ([!])) and the problem reduces to 
studying Brownian motion on the unit sphere, which can 
be handled by standard operator techniques familiar from 
quantum mechanics. 

Let us suppose to begin with that the initial ({a — 
7r| s =o) and final (tg = 4§ \s=l) tangent vectors are held 
fixed. Q{r) has the path integral representation 

AAlp^le- 1 /'^^ 2 ^^^) 2 ^^^^^^) (2) 

where Af is a normalisation constant. Instead of Q(r) we 
focus on the quantity P(z) = J drQ(r)5(r3 — z), which 
is Q(r) integrated over a plane of constant z. Note that 
P(z) and Q(r) vanish when the modulus of their argu- 
ments exceeds L. The generating function of P{z) is de- 
fined as P(f) = J^ L dze fz / L "P{z), where L p = A/k B T. 
Performing the elementary integrations involving 5- func- 
tions we find that P(f) can be expressed as Z(f)/Z(0), 
where Z(f) has the path integral representation 

Z{f)=N J V [i( s )] e ~ Lp/2[ fa {di/ds)2ds] e [f/Lp fo Lds] 

(3) 

Making the change of variable r = s/L p , we arrive at the 
expression 

Z{f) = M J V[t{j)]e~ ^ dr[l/2(dt/drf-fL] ^ 

where (3 = L/L p . Eq. (Q) can be interpreted as the path 
integral representation for the kernel of a quantum par- 
ticle on the surface of a sphere at inverse temperature 
(3. Thus we can express Z(f) as the quantum amplitude 
to go from an initial tangent vector Ia to a final tan- 
gent vector is in imaginary time (3 in the presence of an 
external potential —/cos 0: 

Z(f)=Y,e- [PEn] r n (iA)i>n(tB). (5) 



Here {ip n (t)}, is a complete set of normalized eigenstates 

of the Hamiltonian H = — /cos (9 and E n are the 

corresponding eigenvalues. 

In this paper we focus on the situation where the 
boundary tangent vectors are unconstrained, i. e they are 
integrated over with uniform weight. For free boundary 
conditions the situation is spherically symmetric and Q 
depends only on r = \r\ and not on r and we can write 
Q(r). The probability distribution for the end to end 
distance r is given by S(r) — 47rr 2 Q(r), as can be seen 
by integrating Q(r) over a sphere of radius r. As we 
did above, we can integrate Q{r) over a plane of fixed 
z:P(z) = J dfQ(r)6(r 3 — z). P(z) is the probability dis- 
tribution for the z co-ordinate of the tagged end. Both 
S(r) and P(z) are experimentally accessible quantities 
and they are integrals of the spherically symmetric func- 
tion Q(f) over two dimensional surfaces. Using tomo- 
graphic techniques (reconstruction of a function from a 
knowledge of its integral over two dimensional slices) one 
can deduce jlTj the relation S(r) = —2rdP(r)/dr, where 
P(r) is P{z) with its argument replaced by r. 

For free boundary conditions (|5|) can be written as 
a "vacuum persistence amplitude" 

Z(/)=<0|exp-/3 J ff / |0> (6) 

where Hf = — -^V 2 — fcosO is the Hamiltonian of 
the rigid rotor Q in a potential and |0 > is the 
ground state of the free Hamiltonian Ho = — ^V 2 . 
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FIG. 1. G(f) as a function of / for /3 = L/L p = 1. Also 
plotted are approximate analytic forms valid in the small 
(open circles) and large (filled circles) force regimes. 

By choosing a basis in which Hq is diagonal we find 
that H is a symmetric tridiagonal matrix with diagonal 
elements Hn = 1(1 + 1)/ '2 and superdiagonal elements 
Hi i+i = f(l + l)i/l/((2Z + l)(2/ + 3)). Upto this point 
the treatment is completely analytical. To evaluate Eq. 

we need to use numerical methods. Hf is really an 
infinite matrix, but we truncate it to NXN size, numer- 
ically evaluate it (using Mathematica p8j) and adjust 
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the cutoff N until the answer stabilizes to desired accu- 
racy. From this we deduce all the properties of the model, 
to an accuracy limited only by computational power. 
The form of G(f) = -l//31ogZ(/) is shown in Fig. 1 
along with physically motivated approximate analytical 
forms valid in the small ( G(f) = Ci(/3)/ 2 ) and large 
(G(f) = -f + v7- log//(2/3) + C 2 (/3)) force regimes. 
More terms can be computed, but this already gives a 
fair fit. 

From P(f) (which is equal to Z(f) since Z(0) = 1), 
it is possible to compute P(z) by performing the inverse 
Laplace transform. (Numerically it is more convenient 
to use the inverse Fourier transform by going to imag- 
inary /). The results are shown in Figure 2. For con- 
venience, we set L p = 1 so that (3 = L and plot all 
figures in terms of scaled variables, f = zj '/3, p = rj (3. 
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FIG. 2. The distribution P(Q of scaled extension £ = z//3 
for (3 = L/Lp equal to 1, 3.85 and 10. 
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FIG. 3. The distribution S(p) of scaled end-to-end dis- 
tances p = r//3 for f3 — L/L p equal to 1, 3.85 and 10. Dots 
show simulation data taken from Ref. [9]. 

From the relation S{p) = —2p-^P(p) we compute the 
distribution of end-to-end distance. These are displayed 



in Fig. 3. We have checked that these graphs quan- 
titatively agree (to within the errors of the simulation 
data) with the published plots of fL2fl . Notice that P(() 
and S(p) both have a single maximum and the corre- 
sponding free energies have a single minimum. How- 
ever, for a range of (3 near 3.8, Q(p) develops a dou- 
ble humped form, reflecting the existence of two stable 
free energy minima resulting in a "first order transition" , 
where the quotes signify that this is not a true phase 
transition due to finite size effects. This feature was first 
noticed in computer simulations of the WLC model. 
Our theoretical work confirms the results of simulations 
presented in ||. The form of Q(p) is plotted in Fig.4 
along with the results of computer simulations from []9| . 



0.8 




Q(P) 






0.6 




H 


0.4 


P=3.85 \\ / 




0.2 







0.2 0.4 0.6 0.8 1 1.2 
P 



FIG. 4. The function Q(p) for /3 = L/L p = 1,3.85 and 5. 
Results of a simulation from Ref. [10] are also plotted on the 
curve for j3 = 3.85. 




FIG. 5. The mean extension f as a function of / (thick line) 
for P = L/Lp = 10. Also shown on either side are the root 
mean square fluctuations (thin lines) of the extension about 
its mean value. 

A property of direct experimental interest is the force- 
extension relation (FER) . We work in the constant / en- 
semble and in Fig. 5 plot the scaled mean extension ((f) 
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(defined by ( = -dG(f)/df). 

Since we are dealing with a system of finite size, we 
expect that the extension £ will fluctuate about its mean 
value C- The theoretically expected root mean square 

value of these fluctuations AC = J 1/ ' f3dQ/df of Q is 
shown in Fig. 5. These fluctuations clearly vanish in 
the limit of infinitely long polymers |jj . 

Although the WLC model has been known for around 
fifty years |llj , there does not appear to be a closed form 
analytic solution in terms of elementary functions. This 
paper presents a numerical solution to the WLC, which, 
given the power of modern personal computers is as use- 
ful as an exact analytic form. Using the techniques out- 
lined in this paper one could work out the predictions 
of the WLC model to any desired accuracy, for example, 
experimental accuracy. Our work provides predictions of 
force-extension curves for all lengths which can be tested 
against experiments. All the quantities computed here 
(Q(p),P((),S(p), P(f), G(f) and FER) are experimen- 
tally measurable. Here we briefly go over the pertinent 
experimental realizations of some of these quantities, (a) 
Measurement of Q(p): One can measure Q(p) by attach- 
ing a bead to one end of the molecule and confining the 
bead in a stiff optical trap and by recording the distri- 
bution of location of the other end (tagged with dye) by 
means of a CCD camera, (b) Measurement of P((): P(() 
can be measured by recording the location of the free end 
of the molecule on a given £ plane and focussing all the 
light from the particular £ plane by using a confocal mi- 
croscope. 

In this paper we have used free boundary conditions 
for the tangent vectors. Other boundary conditions can 
also be handled as explained in [jl9| . The boundary con- 
ditions depend on the particular experimental setup. For 
example, if the tangent vectors at the ends are held fixed 
one needs to use S function weights rather than uniform 
ones. Choice of theoretical weights consistent with ex- 
perimental boundary conditions is particularly crucial in 
the context of short polymers. 

The FERs predicted by the model depend on the en- 
semble in which the calculation is done. Depending on 
the experimental situation one should use an ensemble 
in which one of z, r, r or their conjugate forces /, / r , / is 
held constant. As an example we display in Fig. 5 the 
FER in the constant / ensemble. The force extension 
relations in this ensemble are monotonic for all values 
of (3. (In contrast, the FER in the constant p ensemble 
is non-monotonic § in the (3 range where the function 
Q(p) is double humped.) As mentioned earlier, since we 
are dealing with a finite system, which is not near the 
thermodynamic limit, equivalence between conjugate en- 
sembles is not assured |Tc|| . This is due to fluctuations 
around mean values which are not negligible for short 
polymers. Our analysis, being exact, correctly takes into 
account such finite size effects. The pure bend WLC 



model has a single parameter (3 = L /L p and predicts not 
only a force extension relation but also the amount of 
theoretically expected noise on this curve. In an experi- 
ment one can expect to see this noise over and above any 
instrumental noise that may be present in the system. 

Self-avoidance is a feature present in real polymers, 
which has not been taken into account in this analysis. 
Such effects are more important for flexible polymers. 
Self-avoidance is difficult to handle analytically and is 
one of the important directions for future work. We hope 
that the results of this paper will stimulate a detailed and 
quantitative comparison between the predictions of the 
WLC model and experiments, and lead to an improved 
understanding of the elasticity of semi-flexible polymers. 

Note added: After this manuscript was submitted 
for publication, we learned of closely related work by 
Stepanow and Schutz |po[ . 
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